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In  this  work,  using  an  in-house  made  Loschmidt  diffusion  cell,  we  measure  the  effective  coefficient  of 
dry  gas  (O2-N2)  diffusion  in  cathode  catalyst  layers  of  PEM  fuel  cells  at  25  °C  and  1  atmosphere.  The 
thicknesses  of  the  catalyst  layers  under  investigation  are  from  6  to  29  pim.  Each  catalyst  layer  is  deposited 
on  an  A1203  membrane  substrate  by  an  automated  spray  coater.  Diffusion  signal  processing  procedure  is 
developed  to  deduce  the  effective  diffusion  coefficient,  which  is  found  to  be  (1.47  ±  0.05)  x  10-7  m2  s-1 
for  the  catalyst  layers.  Porosity  and  pore  size  distribution  of  the  catalyst  layers  are  also  measured  using 
Hg  porosimetry.  The  diffusion  resistance  of  the  interface  between  the  catalyst  layer  and  the  substrate  is 
found  to  be  negligible.  The  experimental  results  show  that  the  O2-N2  diffusion  in  the  catalyst  layers  is 
dominated  by  the  Knudsen  effect. 

Crown  Copyright  ©  2010  Published  by  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

A  cathode  catalyst  layer  (CCL)  in  a  proton  exchange  membrane 
(PEM)  fuel  cell  is  a  porous  thin  coating  with  the  porosity  of  30-60% 
and  pore  size  distribution  from  several  nanometers  to  100  nm.  The 
porosity  allows  the  transport  of  reactant  oxygen  into  the  CCL  and 
the  transport  of  product  water  vapour  out  of  the  CCL.  Thus  oxygen 
reduction  reaction  can  happen  through  the  whole  depth  of  the  CCL, 
maximizing  current  density.  Oxygen  concentration  in  air  is  about 
20%  and  02  transports  into  the  CCL  mainly  through  a  diffusion  pro¬ 
cess.  Consequently,  the  effective  gas  diffusion  coefficient  (EGDC) 
of  the  CCL  is  a  very  important  parameter,  which  affects  the  uni¬ 
formity  of  oxygen  reduction  reaction  through  the  whole  CCL,  the 
CCL  lifetime,  and  the  power  density  of  PEM  fuel  cells  [1  ].  Therefore, 
it  is  in  great  demand  to  theoretically  calculate  this  parameter  and 
experimentally  measure  it. 
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Theoretically  the  EGDC  of  a  porous  material  is  correlated  to 
the  corresponding  gas  diffusion  coefficient  D  through  empty  space, 
namely  the  bulk  diffusion  coefficient,  and  some  statistical  param¬ 
eters  of  porous  microstructure,  such  as  porosity.  Some  theoretical 
formulas  describing  the  relation  between  the  EGDC  and  porosity 
are  available  [2-8].  Among  the  formulas,  the  Bruggeman’s  formula 
[4]  is  very  frequently  used  for  calculating  the  EGDC  of  a  CCL  But 
the  accuracy  of  the  calculated  EGDC  is  often  inadequate  and  even 
misleading  when  pore  size  distribution  of  the  CCL  falls  in  the  region 
of  less  than  1  micron,  i.e.,  the  Knudsen  effect  must  be  considered. 
In  order  to  quantify  the  Knudsen  effect,  a  3D  microstructure  based 
method  for  calculating  the  EGDC  was  developed,  and  a  formula 
relating  the  EDGC  to  porosity  and  Knudsen  effect  was  proposed  by 
Mu  et  al  [2]. 

Experimental  measurement  techniques,  e.g.,  gas  chromatogra¬ 
phy,  nuclear  magnetic  resonance  (e.g.,  PFG  NMR),  and  diffusion  cell 
methods,  have  been  developed  to  determine  the  EGDCs  of  porous 
materials  [9-13].  Kramer  and  co-workers  measured  the  EGDC 
in  carbon  paper  using  electrochemical  impedance  spectroscopy 
[14,15].  Zhang  et  al.  [9]  used  a  Wicke-Kallenbach  diffusion  cell 
to  measure  the  EGDC  of  a  catalyst  monolith  washcoat.  A  closed- 
tube  method  with  a  Loschmidt  diffusion  cell  is  considered  as  one 
of  the  most  reliable  methods  to  determine  binary  diffusion  coeffi¬ 
cients  of  gases  [16-19].  Using  a  Loschmidt  cell,  bulk  binary  diffusion 
coefficients  of  O2-N2  were  precisely  measured  under  the  experi¬ 
mental  conditions  of  different  temperatures  (25-80  °C)  and  relative 
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humidity  values  (0-80%)  [18];  the  EGDCs  of  gas  diffusion  layers 
of  PEM  fuel  cells  were  also  studied  at  different  temperatures  [11]. 
Moreover,  the  measured  EGDC  of  a  porous  sample  of  a  stainless 
steel  film  with  simple  straight  pores  was  found  to  be  in  good 
agreement  with  the  result  of  numerical  computation  of  three- 
dimensional  mass  diffusion  through  the  sample  [19],  exhibiting  the 
usefulness  of  the  Loschmidt  cell  in  studying  porous  material. 

Although  the  EGDCs  of  the  CCLs  in  PEM  fuel  cells  are  in  great 
demand,  and  measurement  techniques  of  the  EGDC  have  been 
developed,  little  literature  can  be  found  on  the  measured  EGDCs 
of  the  CCLs  of  PEM  fuel  cells.  Over  the  past  several  years,  we  have 
been  trying  to  correlate  the  CCL  microstructure  with  its  EGDC  and 
then  its  performance.  As  a  part  of  this  effort,  we  have  developed 
a  Loschmidt  cell  method  for  measuring  the  EGDCs  of  CCLs.  In  this 
paper,  we  present  the  measurement  method  first  and  then  the  mea¬ 
sured  EGDC  of  eight  CCL  samples,  which  are  different  in  thickness 
from  6  to  29  pum  but  are  fabricated  with  the  same  composition 
and  the  same  deposition  conditions,  allowing  us  to  consider  the 
microstructures  of  these  CCLs  to  be  similar  in  a  statistical  sense. 
Each  sample  is  deposited  on  an  A1203  membrane  substrate.  In  order 
to  measure  the  EGDC  of  eight  CCLs,  the  bulk  diffusion  coefficient 
and  the  EGDC  of  the  substrate  have  to  be  measured  first.  Our  mea¬ 
surement  indicates  that  the  EGDC  is  one  order  of  magnitude  lower 
than  the  Bruggeman  formula  [4]  gives,  meaning  that  the  Knudsen 
effect  is  dominant  and  must  not  be  neglected.  This  result  confirms 
our  previous  calculations  of  the  EGDCs  of  CCLs  [2]. 

2.  Experimental 

2.  I .  Experimental  setup 

The  in-house  made  Loschmidt  diffusion  cell  comprised  two 
chambers,  separated  by  a  valve  (atz  =  0),  as  schematically  shown  in 
Fig.  1.  The  chambers  were  maintained  at  25  °C  and  1  atmosphere 
in  this  work.  At  the  beginning  of  a  measurement,  the  top  cham¬ 
ber  was  filled  with  N2,  while  the  bottom  chamber  was  full  of  02. 
When  the  valve  opened  the  binary  gas  diffusion  started,  and  an 
optic  fibre  oxygen  sensor  located  at  (z=Z)  measured  the  oxygen 
concentration  value  in  the  top  chamber.  The  diffusion  coefficient 
was  deduced  from  the  oxygen  concentration  evolution  with  time 
in  the  top  chamber  detected  by  the  oxygen  sensor.  A  computer 
controlled  the  whole  experimental  procedure,  including  gas  fill¬ 
ing  the  chambers,  signal  collection  from  the  oxygen  sensor,  and 
temperature  maintenance  of  the  cell.  The  detailed  description  of 
the  experimental  setup  and  experimental  procedure  can  be  found 
in  Refs.  [11,18].  In  this  work,  the  ball  valve  in  Refs.  [11,18]  was 
replaced  by  a  flat  one,  and  Z=  1.595  x  10-2  m. 

2.2.  Processing  of  diffusion  signals 

In  this  work  for  the  Loschmidt  diffusion  cell  shown  in  Fig.  1, 
the  02-N2  diffusion  coefficient  can  be  attained  from  the  oxygen 
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Fig.  1.  A  schematic  diagram  of  a  Loschmidt  cell  consisting  of  two  chambers:  bottom 
and  top.  A  flat  valve  was  located  at  z  =  0.  An  optic  fibre  oxygen  sensor  was  poisoned 
at  z  =  Z.  The  diameter  of  the  oxygen  sensor  was  300  |jim. 


concentration  evolution  in  the  top  chamber  based  on  the  one¬ 
dimensional  Fields  second  law  of  diffusion  [16,17,20]: 

dC(z,  t)  _Dd2C(z,  t) 


with  the  initial  conditions  and  the  impermeable  boundary  condi¬ 
tions  atz  =  ±L/2: 


c  =  cb 

(-h2<  »• '-») 

|,  (2) 

c  =  ct 

(°<zs+t,  '-») 

1,  (3) 

(  ar\ 

=  0  (t>0). 

z=±L/2 

(4) 

C(z,t)  is  the  oxygen  concentration,  D  the  diffusion  coefficient.  L  is 
the  length  of  the  Loschmidt  cell.  Cb  and  Ct  are  the  concentrations  of 
02  at  t  =  0  in  bottom  and  top  chambers,  respectively.  The  solution 
of  Eq.  (1)  is  [17-19] 


C  (z,  t)  = 


1 


Cb  +  Q-(Q-Ct)^er/( 


2  VDt 


(t<0.1r).(5) 


er/(x)  is  the  error  function,  and  r  =  L2/  (nD)  is  the  characteristic  dif¬ 
fusion  time.  Refs.  [17,18]  exhibit  that  for  a  short  period  t  <  0.1  r,  Eq. 
(5)  well  describes  the  bulk  diffusion  (diffusion  without  any  porous 
sample  in  the  diffusion  cell)  of  a  binary  gas  mixture.  The  bulk  dif¬ 
fusion  coefficient  D  can  be  deduced  by  least-square  curve  fitting 
Eq.  (5)  to  the  curve  of  oxygen  concentration  evolution  in  the  top 
chamber  sensed  by  the  oxygen  sensor. 

In  the  presence  of  a  porous  sample  located  between  z  =  0  and 
z=Z  in  the  top  chamber  in  Fig.  1,  an  equivalent  resistance  Req  to 
binary  gas  diffusion  may  be  introduced  [9]. 


^  ~  Wa. A  ’ 


(6) 


where  Deq  is  the  equivalent  diffusion  coefficient  of  the  medium 
from  the  z  =  0  to  z=Z,  and  A  is  the  cross-sectional  area  available  to 
diffusion.  In  this  case,  the  D  in  Eq.  (5)  should  be  replaced  by  Deq.  The 
oxygen  concentration  evolution  curve  after  diffusion  through  the 
porous  sample  can  be  measured  with  the  oxygen  sensor,  and  the  Deq 
can  be  deduced  with  Eq.  (5)  by  the  foregoing  curve-fitting  method. 
The  concept  of  the  equivalent  resistance  originates  from  electrical 
conduction,  and  sometimes  heat  conduction  analysis  adapts  the 
concept.  In  this  method,  the  object  is  considered  to  be  comprised  of 
a  number  of  individual  resistances,  which  can  be  connected  to  form 
an  overall  circuit  [9].  In  the  presence  of  a  porous  sample  located 
between  z  =  0  and  z = Z,  the  equivalent  resistance  to  gas  diffusion  is 
formed  by  bulk  diffusion  resistance  and  porous  sample  resistance, 
and  they  are  connected  in  series. 

In  the  case  of  a  single  layer  porous  sample  of  thickness  l  and 
effective  diffusion  coefficient  Dejj,  the  equivalent  resistance  can  be 
considered  the  sum  of  the  bulk  diffusion  resistance  and  the  diffu¬ 
sion  resistance  through  the  porous  sample. 


_  Z  _Z-l  l 
eq  ~  ~D^A  ~  ~DA  +  De(fA ' 


(7) 


The  Deff  can  be  calculated  with 


°eff  Z/Deq  -  (Z  -  l)/D ' 


(8) 


If  the  sample  is  of  double  layers,  e.g.,  a  CCL  deposited  on  a  substrate, 
the  equivalent  resistance  Rceq  may  be  expressed  as 


Rceq  — 


z 

DceqA 


Z  —  Is  —  lc  Is  lc 
DA  +  DseffA  +  DceffA 


(9) 
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Here  ls  and  lc  are  the  thickness  of  the  substrate  and  catalyst  layer, 
respectively.  Dseff  and  Dceff  are  the  effective  diffusion  coefficients  of 
the  substrate  and  catalyst  layer,  respectively.  Dceq  is  the  equivalent 
diffusion  coefficient  in  the  presence  of  the  CCL  and  the  substrate. 
The  effective  diffusion  coefficient  of  the  catalyst  layer  Dceff  can  be 
found  with  Eq.  (10): 

Dceff  =  ZjlTTq  -  (Z  -  /s  -  M/D  -  UjD^f  ’  (10) 

with  Eqs.  (5),  (8)  and  (10),  the  effective  diffusion  coefficient  of  a 
porous  sample,  e.g.,  a  CCL,  can  be  determined. 

2.3.  CCL  sample  preparation 

The  CCL  samples  contained  30  wt%  Nation  ionomer  (5  wt%  solu¬ 
tion,  EW1000,  Alfa  Aeser).  The  Pt/C  catalyst  used  in  this  work  was 
46  wt%  Pt/graphitized  carbon  (Tanaka  Kikinzoku  Kogyo).  The  base¬ 
line  ink  formulation  contains  1:1  mixture  of  methanol/water.  The 
mixture  was  sonicated  in  a  water  bath  at  room  temperature  for 
30  min.  Nation  solution  was  added  drop-wise  to  the  mixture.  The 
suspension  was  homogenized  for  60  minutes.  The  ink  was  ready 
to  be  used  for  deposition.  An  automated  spray  coater  (EFD-Ultra 
TT  series)  was  employed  to  deposit  CCLs  onto  the  A1203  membrane 
substrates,  Anopore  Inorganic  Membranes  (Anodise  25,  Whatman). 
Various  Pt  loadings,  0.2,  0.3,  0.4,  0.6  and  0.8  mg  cm-2,  correspond¬ 
ing  to  about  6-29  p,m  in  thickness,  were  obtained  by  changing  the 
number  of  deposition  passes.  The  total  of  eight  CCL  samples  with 
different  thicknesses  were  investigated  in  this  work.  The  porosity 
and  pore  size  distribution  of  the  catalyst  layer  were  obtained  by  Hg 
porosimetry  using  the  same  procedure  as  previously  reported  [21  ]. 

3.  Results  and  discussion 

3.1.  Substrate 

In  this  work,  each  CCL  sample  under  investigation  was  deposited 
in  a  piece  of  substrate.  The  functional  requirements  of  the  substrate 
are:  (1)  possessing  well-defined  pore  structure  with  no/minimum 
lateral  diffusion  within  the  substrate,  (2)  exhibiting  considerably 
higher  diffusion  coefficient  than  that  of  the  catalyst  layer  such  that 
the  limiting  transportation  is  not  arising  from  the  substrate,  and 
(3)  having  pore  diameters  small  enough  such  that  the  catalyst  layer 
deposited  would  not  penetrate  into  the  substrate,  and  thus  a  well- 
defined  catalyst  layer  thickness  can  be  determined.  In  addition,  the 
substrate  must  be  chemically  inert  to  common  solvents  used  during 
a  catalyst  layer  deposition  (e.g.,  low  b.p.  alcohols  and  water). 

We  found  that  AI2O3  membrane  was  a  good  candidate  of 
the  substrate  satisfying  the  requirements.  The  membrane  is 
hydrophilic  and  compatible  with  most  solvents  and  aqueous  mate¬ 
rial.  According  to  the  vendor,  this  material  has  a  non-deformable 
honeycomb  pore  structure  with  no  lateral  crossovers  between  indi¬ 
vidual  pores.  Fig.  2  is  an  SEM  image  of  the  cross  section  of  the 
membrane  “Anodise  25”  from  Whatman  International  Ltd.  The 
image  shows  how  the  pore  size  changes  from  circa  200  nm  on  the 
bottom  surface  to  about  20  nm  on  the  top  surface.  Table  1  lists  the 
geometrical  parameters  of  the  Anodise  25  quoted  from  the  vendor 
and  measured  in  this  work. 

Three  Anodise  25  substrates  were  measured,  and  each  substrate 
was  measured  more  than  five  times.  Fig.  3  shows  a  typical  curve 
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Fig.  2.  Cross-section  SEM  image  of  the  inorganic  (AI3O2)  membrane  “Anodise  25”. 


of  oxygen  concentration  evolution  in  the  presence  of  an  Anodise 
25  substrate.  Fitting  Eq.  (5)  to  the  concentration  curves,  we  could 
attain  the  equivalent  diffusion  coefficient  Dseq ,  and  the  averaged 
Dseq  of  the  three  substrates  was  Dseq  =  (1.90  ±0.03)  x  10-5  m2  s-1. 
In  the  experiments,  the  exact  time  t0  of  the  beginning  of  the  gas 
diffusion  was  hard  to  determine.  To  find  t0,  the  t  in  Eq.  (5)  was 
replaced  by  t  - 10  when  curve  fittings  were  performed.  With  Eq. 
(8),  the  EGDC  Dseff  of  each  substrate  was  calculated,  and  the  arith¬ 
metic  average  of  the  EGDCs  of  the  three  substrates  was  found  to  be 
Dseff=(  1.3  ±0.3)  x  10-6  m2  s-1.  In  the  calculation  using  Eq.  (8),  bulk 
02-N2  diffusion  coefficient  D  =  2.02  x  10-5  m2  s-1  [18]  was  used. 

3.2.  Pore  size  distribution  of  the  CCLs  and  porosity 

The  peak  pore  size  and  pore  size  distribution  (PSD)  of  the  CCLs 
were  examined  by  Hg  porosimetry,  as  shown  in  Fig.  4.  In  general, 
the  pore  size  distribution  of  these  CCLs  ranges  from  10  to  200  nm, 
which  includes  meso  (2-50  nm  diameter)  and  macro  (>50nm) 
pores  as  defined  by  the  International  Union  of  Pure  and  Applied 
Chemistry  (IUPAC).  The  peak  pore  size  for  these  CCLs  is  about  50  nm. 
The  cumulative  intrusion  curve  represents  the  integral  volume  of 


Fig.  3.  The  oxygen  concentration  evolution  in  the  top  chamber,  in  the  presence  of 
an  Anodise  substrate,  and  a  CCL  (29  p,m)  on  the  top  of  the  substrate,  respectively. 


Table  1 

The  geometrical  parameters  of  Anodise  25  inorganic  (A1203)  membranes. 


Average  thickness 

Diameter  with  ring 

Membrane  diameter 

Pore  size 

Porosity 

From  the  vendor  60  p.m 

Measured  with  mercury  porosimetry  in  this  work 

25  mm 

21  mm 

0.02  |jim 

25-50% 

28% 
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Fig.  4.  Typical  pore  size  distribution  of  the  CCLs. 


Fig.  5.  Numerical  simulation  of  dDceffldDceq  using  Eq.  (10). 


Hg  penetrated  into  the  pores  with  diameters  ranging  from  1  |±m 
to  5  nm.  The  porosities  of  the  CCLs  with  different  thickness  were 
subsequently  calculated  to  be  around  30-40%. 

3.3.  The  effective  gas  diffusion  coefficient  of  the  CCL 

The  thickness  of  CCL  samples  varied  from  circa  6  p,m  to  about 
29  |±m,  and  eight  samples  of  different  thicknesses  were  studied  in 
this  work,  as  shown  in  Table  2.  Each  sample  was  measured  more 
than  five  times.  A  typical  curve  of  oxygen  concentration  evolution 
in  the  presence  of  a  CCL  and  the  substrate  was  shown  in  Fig.  3. 
After  the  determination  of  the  equivalent  diffusion  coefficient  using 
Eq.  (5),  the  EGDC  of  each  sample  could  be  deduced  with  Eq.  (10). 
Table  2  exhibits  that  the  DceJgr  does  not  change  significantly  with 
the  thickness  of  CCLs,  and  the  arithmetic  average  EGDC  of  the  eight 
CCLs  was  found  to  be  DceJ^=(1.5±0.1)  x  10-7  m2  s-1. 

It  was  observed  that  the  EGDC  of  thin  (<10pJm)  CCL  samples 
is  more  sensitive  to  the  fluctuation  of  Dceq  than  that  of  the  thick 
(>10  pan)  CCL  samples,  resulting  in  a  larger  deviation  of  the  EGDC 
of  the  thin  CCL  samples.  From  Eq.  (10),  one  may  have 

ADceff  =  ^p^ADceq.  (11) 

UL>ceq 

Fig.  5  reveals  that  dDceffldDceq  increases  when  the  CCL  thickness 
decreases,  giving  the  explanation  of  the  observed  phenomenon. 

To  reduce  the  influence  of  the  fluctuation  on  the  EGDC,  the  aver¬ 
aged  EGDC  could  be  determined  by  a  least-square  curve  fitting  of 
the  following  equation,  derived  from  Eq.  (9), 

Z 

^  (Z  —  ls  —  lc  )/D  +  ls/ Dsejj  +  lc/ Dceff 

to  the  experimental  data  presented  in  Fig.  6.  With  the  correla¬ 
tion  coefficient  R  =  0.984  of  the  curve  fitting,  the  EGDC  thus  attained 
was  Dceff=  (1.47  ±0.05)  x  10-7  m2  s-1,  consistent  with  the  forego¬ 
ing  arithmetic  average  DceJgr=(1.5 ±0.1)  x  10-7  m2  s-1.  Moreover, 


Table  2 

Measured  equivalent  and  effective  binary  gas  diffusion  coefficients  of  catalyst  layers. 


lc(  ixm) 

Dceq(  10-5  m2  S'1) 

Dceff  i  10“7  m2  s-1) 

6 

1.81  ±0.03 

1.36 

9 

1.81  ±0.04 

1.67 

10 

1.73  ±0.04 

1.24 

12 

1.75  ±0.02 

1.50 

14 

1.72  ±0.02 

1.47 

20 

1.67  ±0.03 

1.62 

23 

1.62  ±0.03 

1.50 

29 

1.54  ±0.01 

1.43 

the  EGDC  obtained  using  the  curve  fitting  method  has  one  more 
significant  digit  than  the  arithmetic  average,  improving  the  mea¬ 
surement  precision. 

In  Eqs.  (9),  (10)  and  (12)  the  possible  diffusion  resistance  of 
the  interface  between  the  substrate  and  the  CCL  on  the  top  of  the 
substrate  is  not  taken  into  account.  To  find  the  possible  diffusion 
resistance,  we  rewrite  Eq.  (9)  as 


Z 

Dceq 


Z-ls-lc  |  *s  |  lc 
D  DSeff  Dceff 


lc 


(13) 


Eq.  (13)  exhibits  that  RceqA  is  a  linear  function  of  the  CCL  thickness 
Zc;  the  slope  of  RceqA  is  basically  decided  by  Dceff,  and  the  intercept 
should  be  Z/Dseq  if  there  is  no  diffusion  resistance  at  the  interface. 
Fig.  7  presents  a  least-square  linear  fitting  to  the  expermental  data 
of  ZIDceq.  For  a  linear  fitting  there  exists  a  best  approximation  [22]. 
The  intercept  of  the  linear  fitting  was  found  to  be  8.33  x  103  m-1  s. 
With  Dseq  =  (1.90  ±0.03)  x  10-5  m2  s-1,  Z\Dseq  is  in  the  range  of 
8.26-8.53  x  103  nrr1  s.  The  intercept  found  is  within  the  range,  indi¬ 
cating  that  the  diffusion  resistance  of  the  interface  is  negligible. 


Fig.  6.  Experimental  data  of  the  equivalent  diffusion  coefficients  of  the  eight  CCLs 
of  different  thicknesses  and  the  least-square  curve  fitting  Eq.  (12)  to  the  experi¬ 
mental  data:  Z=  1.595  x  10-2  m,DseJ^=(1.3  ±0.3)  x  10-6  m2  s-1,and  /s  =  6.0  x  10_5m 
were  used  in  the  fitting.  The  correlation  coefficient  between  the  fitting  and  the 
experimental  data  R  =  0.984  was  found. 
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Fig.  7.  A  best  linear  fitting  to  the  expermental  data  of  RceqA  (i.e.,  Z/Dceq). 


3.4.  The  Knudsen  effect 

The  Bridgeman  formula,  Dceff=Ds1-5,  has  been  widely  used 
for  calculating  the  EGDC  of  catalyst  layers  of  PEM  fuel  cells. 
Here  s  is  porosity.  Through  these  experimental  measure¬ 
ments  of  the  8  CCL  samples,  we  found  that  Bruggeman’s 
formula  overestimates  the  EGDC  of  the  CCL  samples.  The 
porosity  of  the  eight  samples  ranges  from  30  to  40%,  so 
the  average  EGDC  of  the  8  samples,  calculated  based  on 
Dce//=D£15,  is  about  Dceff=  2.02  x  10-5(0.35)15  =4.18  x  10"6  m2 
s-1.  Comparing  with  the  measured  EGDC  of  the  CCL  samples, 
Dce//=(1.5  ±0.1)  x  10-7  m2  s-1,  the  EGDC  of  the  CCLs  from  Brugge¬ 
man’s  formula  is  one  order  higher  than  the  measured  one.  This 
considerable  difference  is  due  to  the  fact  that  Bruggeman’s  formula, 
Dceff =De1-5,  fails  to  take  into  account  the  Knudsen  effect.  Generally 
speaking,  the  Knudsen  effect  is  not  negligible  when  pore  size  is  as 
small  as  one  micrometer,  and  becomes  dominant  when  pore  size 
is  less  than  100  nm.  As  can  be  seen  from  Fig.  4,  the  pore  size  distri¬ 
bution  of  the  CCL  samples  ranges  from  10  to  200  nm  with  the  peak 
pore  size  around  50  nm,  meaning  that  the  Knudsen  effect  has  to  be 
considered. 

Mu  et  al.  [2]  proposed  a  method  for  calculating  the  EGDC  with 
the  Knudsen  effect  included  using  a  finite  element  method.  In  Mu’s 
method  [2],  the  entire  pore  space  in  a  porous  layer  is  discretized 
into  a  number  of  small  tubes  with  different  diameters.  The  porosity, 
the  average  diameter  and  the  standard  deviation  of  diameters  of  the 
small  tubes  can  be  specified.  In  a  single  cylindrical  pore,  the  effec¬ 
tive  diffusion  coefficient,  Dejj,  is  represented  by,  1  /Dejj=  1  IDb  + 1  /Dk, 
where  Db  is  the  bulk  diffusion  coefficient,  and  Dk  is  the  Knudsen 
diffusion  coefficient,  which  is  given  by  Dk  =  48.5dp^/T/M.  Here  dp 
is  pore  diameter,  T  is  temperature  and  M  is  molecular  weight.  Mu’s 
results  show  that  when  the  Knudsen  effect  is  counted,  the  effective 
gas  diffusion  coefficient  of  the  catalyst  layer  with  mean  pore  size 
50  nm  and  porosity  40%,  is  about  15%  of  the  effective  gas  diffusion 
coefficient  without  the  Knudsen  effect  included.  This  is  basically  in 


agreement  with  the  effective  gas  diffusion  coefficients  measured 
experimentally  in  this  paper. 

4.  Conclusions 

An  in-house  made  Loschmidt  diffusion  cell  was  employeed  to 
measure  the  effective  coefficient  of  binary  gas  diffusion  through 
cathode  catalyst  layers  of  PEM  fuel  cells.  Each  CCL  under  investi¬ 
gation  was  deposited  on  a  substrate  of  A1203  membrane,  and  the 
thicknesses  of  the  eight  CCLs  varied  from  6  to  29  p,m.  With  the  diffu¬ 
sion  signal  processing  procedure  developed  in  this  work,  we  found 
that  the  EGDC  of  the  substrate  was  one  order  higher  than  that  of 
the  CCLs,  exhibiting  that  the  binary  gas  diffusion  through  the  CCL- 
substrate  assembly  was  dominated  by  the  CCL.  Also,  the  diffusion 
resistance  of  the  interface  between  the  CCL  and  the  substrate  did 
not  play  a  significant  role  in  the  diffusion.  The  EGDC  of  the  CCLs  was 
found  to  be  (1.47  ±0.05)  x  10-7  m2  s-1,  which  is  much  smaller  than 
that  predicted  by  Bruggeman’s  formula.  Our  analysis  shows  that  the 
binary  gas  diffusion  through  the  CCLs  is  dominated  by  the  Knudsen 
effect. 
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